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Abstract 

We study an abstract damped wave equation. We prove that the solution of the damped wave 
equation becomes closer to the solution of a heat type equation as time tend to infinity. As an ap¬ 
plication of our approach, we also study the asymptotic behavior of the damped wave equation in 
Euclidean space under the geometric control condition. 
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1 Introduction 

We study the following abstract damped wave equation 

(1.1) {df+A + Bdt)u = Q. 

Here 9* = ^ and A, B are self-adjoint operators on a Hilbert space H. We assume A is a densely-defined 
positive operator i.e. A > 0. i? is a strictly positive bounded operator on H i.e. B > c Id for a constant 
c > 0. We also assume and B~^/^ are continuous operators on D{A), i.e. there exists a constant 
C > 0 such that for any u S D{A), the following bounds hold 

\\AB^^Mh<C{\\u\\h + \\Au\\h), 

\\AB-^/Mh<C{\\u\\h + \\Au\\h). 

An example of this equation is the usual damped wave equation 

(1.2) {df — A + a{x)dt)u = 0 in R x 

where A is the constant coefficient Laplace operator on R*^, a{x) is a strictly positive definite smooth 
function with bounded derivatives. It is known that eventually the solution of HH) are close to the 
solution of a heat type equation. Here we see this by heuristic argument. The stationary equation 
is closely related to the asymptotic profile. So we take u = e^^uo{x) in (11.211 . we have the following 
stationary damped wave equation 

(1.3) (A^ -I- a{x)X — A)uo = 0. 

To the solution of (EH), the influence of large frequency part becomes smaller as time tend to infinity by 
damping. So for small A, if we neglect A^ term in (O, we obtain 

(a(a;)A — A)uo = 0. 

This equation is a stationary equation to the following heat equation 

{a{x)dt - A)uo = 0. 

So we may expect the asymptotic behavior of the solution to EH) is close to the following abstract heat 
equation 

(1.4) B^^^{dt + A)B^/^v = 0. 
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Here A = B The solution of this equation can be written as follows 


Thus we can expect that the solution of the damped wave equation approach to the above solution as 
time tend to infinity. The following statement gives a justification of this argument. 


Theorem 1.1. Letu be the solution to the Cauchy problem of the following abstract damped wave equation 

f {df + A + Bdt)u = 0, 

\ u\t=o = uo,dtu\t=o = ui. 


Then there exists C > 0 such that for any uq S Dip/A), ui G H and t > 1, the following asymptotic 
profile holds 


||u(t) - < Ct-W\uo\\H + WVAuoWh + ||ui||ff). 


Remark 1.2. ITe can get more sharp estimate, see ii’emarfc 15.51 On the other hand, this decay rate is 
optimal, see Proposition 15.61 which is a generalization of the result of m 

In many situation, we impose some assumptions to the initial data. In the case, Theorem II.II does 
not give sufficient information. For example, we consider the constant coefficient damped wave equation 
on 

(Of - A + dt)u = 0. 

The correspondence heat propagator satisfies 

:L^ ^ 


From this decay estimate, we can not get the top term from Theorem 1 1.1 1 if the initial date uo,ui are in 
and d > 3. So we need to impose some more assumption to the equation. For the purpose, we specify 
H = Lf{Tl,p) where (fl,/r) is a cr-finite measure space and assume that {e“‘"^}t>o is also a bounded 
semi-group on 


( 1 . 6 ) 


||e-‘^u|Ui =O(l)||u|Uyt>0. 


We also assume that there exists m > 0 such that the following diffusion estimate holds 

(1.7) ||e-‘^u||i2=0(t—)||u|Ui, t>0. 

Furthermore B^/'^ is also a bounded operator on satisfying 

( 1 . 8 ) ^\\u\W<\\B^/Ml^<C\\u\U^ 

for a constant C > 0. Our model problem satisfies this assumption. The next result is a diffusion 

type estimate for such operators. 

Theorem 1.3. Under the assumptions of Theorem, II. 11 we also assume (11.611 . (11.71) and (11.81) . Then 
there exists a constant C > 0 such that for any uq € D(\/A) D L^, ui € D L^, t > 1, the solution u{t) 
o/dElD satisfies the following asymptotic profile 

\\u{t) - < Ct~^~^{\\uo\\L'^nL^ + II VHmo||l 2 -h ||ui||L 2 nii). 


In the final section, as an application of the above theorems, we consider an example of damped wave 
equation on Euclidean space and we give an asymptotic profile of the solution for a variable coefficient 
damped wave equation, see Theorem 17. II We also treat a perturbation of this case for which the strictly 
positivity of the damping term may not be satisfied but the geometric control condition holds. In this 
case, we can give the similar asymptotic profile, see Theorem 17.91 

Now we remark some related result. On bounded domains or manifolds, under the geometric control 
condition, the exponential energy decay of the solution was proved in p, m- On the other hand, for 
unbounded domains the energy does not decay exponentially. In m, the fact that the decay late of the 
solution is similar to the heat equation is proved to study the semi-linear damped wave equation. The 
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solution of the damped wave equation tend to the solution of the heat equation is well-known for example 
M- This type results may be called diffusion phenomena for the damped wave equation. For variable 
coefficient damping, Y, Wakasugi f[18)i proved the diffusion phenomena for spacial decreasing damping. 
For abstract setting, the diffusion phenomena may be proved in m but at present, the manuscript does 
not published and their result seems to slightly different from ours. 

If the damping term is too weak, for example short range case, the asymptotic behavior of the solution 
is quite different. In the case, the solution tend to the solution of the wave equation and the local energy 
decay is an important problem, see [3]. 


2 Abstract damped wave equation 

We recall the setting of the problem. In this section, we only assume that B is a bounded and non-negative 
operator. We reduce CH) to a first order system. We introduce the following operator 

A = : Ho no, D{A) = {‘(w,u) e Hix H-Aue H,v e D{VA)}. 

Here Bo = Hi x H is the energy space. Hi is the Hilbert space completion of D{A) by the norm 

\\u\\H. = \\VAurH. 

H~^ denotes the dual space of Hi. We note D{\/~A) is naturally embedded in H^. So in what follows, 
we consider D{\/~A) is a subspace oi H^. By using this operator, we can rewrite the wave equation m 
as follows, 



Since we are considering the damped wave equation, we have the following proposition. 

Proposition 2.1. A is a m-dissipative operator. 

Proof. From the definition, we easily prove A is closed. For {u,v) G D{A), we compute as follows 
Re(^‘('u, v))-Ho = A^^^u)h — {Au + Bv, v)h} = —{Bv, v)h < 0. 

Thus A is dissipative. □ 

From this proposition, we have the following resolvent estimate for ReA > 0 

||(A - 

So by Hille-Yosida Theorem, A generates a contraction semi-group 

U{t) = e*-^ : Bo ^ Bo- 

To an initial date *{uo,ui) G Bo, we write the solution as follows 

\u{t),v{t)) = U{tY(uo,ui). 

We restrict the initial date since the energy space Bo is a exotic space. For ‘(mq, ui) S D{A) x D{VA) C 
D{A), since U{ty{uo,ui) is C^, we have the following identity in H^ 

u(t) =uo+ v{s)ds. 

Jo 

By the density argument, the above identity also holds for *(uo,ui) G D(\/A) x H. From the right hand 
side of the above identity.,we can regard u{t) G H. So we introduce the space B = D{i/A) x B by the 
norm 

Uu,v)\\n = {\\u\\l + \\uf^ + \\v\\lV^y 

Since U{t) is a contraction operator, from the above identity, we have the following bound for the 
propagator 

(2.2) \\U{t)\\n^n<C{l + t). 

In this paper, we call u{t) G D{\/~A) for = U{ty{uo,ui), {uo,ui) GH as the solution to (11.511 . 
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3 Resolvent estimate 


We prove some related estimates to the resolvent of A. Probably these estimates are well-known, c.f. [S], 
but for completeness, we give the proof. First we study the following operator which is almost similar to 
the resolvent 

(3.1) E(\) = + BA + A)-\ 

From the positivity of B, we have the following estimate. 

Lemma 3.1. There exists (5 > 0 such that for ReA > —6 and ImA 7 ^ 0, R{A) : H —^ D[A) exists and the 
following bound holds 

(3.2) ||i?(A)||77^77 =0(|ImA|-i). 

Proof. For u G D{A), we have 

(3.3) ((A^ + BA + A)u, u)h = ((|ReAp — |ImAp -1- BReA + A)u, u)h + iImA((2ReA -I- B)u, u)h- 

Taking imaginary part of the above identity and using B > cld we have the following bound 

|((A^ -I- BA + A)u, u)h\ > |ImA((2ReA -I- B)u, u)h\ 

> |ImA|(c-I-2ReA)(zi, w)ij 

for ReA > — c/2. So (A^ -I- BA + A) is injective if ReA is sufficiently small and ImA 7 ^ 0. We can also 
prove (A^ -I- BA + A) is surjective by considering adjoint operator (A^ -I- BA + A) and applying similar 
argument. So the inverse exists and we have the bound. □ 

Taking real part of the identity (13.31) and applying similar argument, we can also prove the following 
lemma. 

Lemma 3.2. //ReA ^ ih/GT}/ Gxzsis (xnjd thG J'oHo'iv'itiq bouTtd holds 

(3.4) ||B(A)||ff^ff = 0(|ReAr'). 

Remark 3.3. From the identity (A^ + BA + A)* = A^ -I- BA + A, we have //(A)* = R{A). So if R{A) 
exists then R{A) exists. 

Next we estimate i?(A) as an operator between different Hilbert spaces. 

Lemma 3.4. ////(A) exists, then the following bound holds 

(3.5) ||/?(A)|| 77 ^ 77 , =0(||//(A)|| 

Proof. For u G H, we have 

\\R{A)u\\%^ = {^/AR{A)u,^/AR{A)u)H 
= {AR{A)u, R{A)u)h 

= {{A + A"^ + BA)R{A)u, RiA)u)H - {{A"^ + BA)R{A)u, R{A)u)h 
= (u, R{A)u)H - ((A^ -h BA)R{A)u, R{A)u )h 
= 0{\\R{A)\\h^h + (|A| + |Ani|B(A)f^^^)hf^. 


□ 

Lemma 3.5. If R{A) exists, then R{A)A : D{A) —> H can be extend to Hi —^ H and the following bound 
holds 

(3.6) \\R{A)A\\h,^h = 0{\\R{A))\\]i^_H + (|A|i/2 + |A|)||//(A)||77^7/). 
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Proof. By R(X)* = R{X) and Lemma 15^ RW is a bounded operator from H to Hi. Since R{X)* = R{X), 
we can regard R{X) as an operator from H~^ to H by duality argument. From Lemma 13.41 we have the 
following estimate 


RiX) : H-^ ^ H = Oi\\RiX)\\]ilH 

We extend A to the operator form Hi to H~^ by using the following identity 

{Au,v)h-^,Hi = i'/Au,'/Av)H- 


From the definition, A satisfies 

A: Hi-^ H-^ =0{1). 

Thus we obtain the lemma. □ 

Lemma 3.6. If R{X) exists for a A yf 0 then R{X){X + B) : D[A) —>• Hi can he extend to Hi —^ Hi and 
the following bound holds 

(3.7) ||i?(A)(A + + l)\\RiX)\\]i!,^ + (1 + |A|)||i?(A)||^^^ + lApi). 

Proof. We will estimate i?(A)(A + B) using the following identity 

R{X){X + B) = \-\r{X)A. 

A A 

For / G D{A), we consider the following equation 

R{X)Af = u. 


This can be written as follows 

Af = (A^ + BX + A)u. 

So we have 

{Au, u)h = (Af, u)h - ((A^ + BX)u, u)h- 
Taking real part, for sufficiently small e > 0, we obtain 

llulll^, <C((|Ap + |A|)||n|||, + ||/||l,J + s||u||l,^. 

Thus 

||i?(A)Al/||^, < cmT + |A|)||i?(A)Al/||^ + ll/II^J. 

Applying (13.611 . the following bound holds 

||i?(A)A/||^, < C((|Ar/2 + |A|)||i?(A)||^/!,^ + (|A| + |Ani|i?(A)||^^^) + l)||/||ffj. 


By density, we conclude the extension and the estimate. □ 

By these lemmas, we can give the existence of the resolvent. 

Proposition 3.7. If R{X) exists for a X ^ 0, then the resolvent (A — A)~^ : Tto —^ Ro exists and the 
following identity holds 


(3.8) 


(A-A)-i 


fR(X)(X + B) R{X)\ 
-R{X)A R{X)X) ' 


The following estimate also holds 

(3.9) ||(A - Ar^Wuo^Uo = 0{{l + \X\-^/^)\\RiX)\\]il^ + (1 + |A|)||i?(A)||^^^) + |A|-i). 


Proof. From Lemma [3.41 Lemma [331 and Lemma 1331 the right hand side of (13.811 can be extend as an 
bounded operator on Ho. By the direct computation, if *{u,v) G D{A) x D{A) then 

/i?(A)(A + B) i?(A) Va _ M ^ \ fu\ _ fu\ 

\ -RiX)A RiX)Xj^^ V R{X)XJ \A X + Bj {vj ~ \vj 
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So by density it is the left inverse. On the other hand for *{u, v) € D{A) x D{-\/~A) 

(\ -Id\ (R{\){\ + B) R{\)\(u\_( Id (A 

A + fij -R{X)A R{X)X) yv) ~ yAR{X){X + B)-{X + B)R{X)A Id) A) 

We use the following trick, for a A ^ 0 

X{AR{X){X + B)-{X + B)R{X)A)u = {AR{X){X^ +BX^A)- (A^ + BX + A)R{X)A)u = 0. 

Thus the density argument, it is also the right inverse. From Lemma |3.41 Lemma 13.51 and Lemma 13.61 
we have the bound. □ 

Remark 3.8. From this proposition, we have the following identities 

((A-y4)”^ *(0 ,u),‘(0,u))ho = X{R{X)u,v)h, 

((A-y4)”^ *(0 ,u),*(z;,0))ho = {R{X)u,v)h^- 

The left hand side is complex analytic since the resolvent is analytic. So the right had side is complex 
analytic so R{X) : H ^ H,H ^ Hi are complex analytic family of operators. 

Applying Lemma HTTl directly to i?(A)(A + B) : H ^ Hi, we also have the following estimates to the 
resolvent. 

Proposition 3.9. If R{X) exists for a X 0, then we have the following estimate 

(3.10) ||(A - A)-y\n^no = (^((|A| + + |Ani|i?(A)||^^^). 

This proposition says if we restrict the domain of the resolvent to R then the singularity, as A tend to 
0, become weaker one. This estimate is important for energy decay. We also have the following estimate. 

Proposition 3.10. If R{X) exists for a X ^ 0, then we have the following estimate 

(3.11) ||(A - A)-y\n^n = 0((1 + |A|-i/^)||i?(A)||^/!,^ + (1 + |A|)||i?(A)||^^^ + |A|-i). 

By Lemma lOI and Lemma [221 we conclude the following estimates of the resolvent. 

Proposition 3.11. If |ReA| < c for sufficient small c > 0 and ImA 7 ^ 0, then the resolvent satisfies 
\\iX-A)-y\no^no=Oi\lmX\-^ + l), 

||(A - Aryin^no = 0{{1 + |A|)|ImA|-i/2 + (|A|i/2 + |An|ImA|-i), 

||(A-Al)-iw^w = 0 (|ImA|-i + l). 

If ReA > 0, the following estimate also holds 

\\iX-A)-y\no^no^Oi\X\-^ + l), 
ll(A-Al) = C>(1 + |A| + |A| ^/^), 

||(A-AI)-1||h^« = (!1(|A|-i + 1). 

4 Abstract energy decay 

We prove the following abstract energy decay estimate to the damped wave equation. The theorem may 
be well known but we prove here for later argument. 

Theorem 4.1. There exists a constant C such that for any *{uo,ui) G R, and t > 1, the following bound 
holds 

\\U{tY{uo,ui)\\no<Cit-^Anuo,ui)\\n). 

Remark 4.2. We give one example which may be useful to understand the proof of this theorem. We 
consider the constant coefficient damped wave equation on R'^ 

{df - A + dt)u = 0. 
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The Fourier transform of the solution is the following form 


u+{f)exp 




+ u_(5) exp 




For simplicity, we assume u+ and U- are in , then U- part decays exponentially so it is no interest. 
If f = 0, it+ part does not decay therefore uniform decay does not occur if we only assume £ 

However if we assume |5|u+(^) £ L^, we have the uniform decay estimate, by using the following 
decomposition. 




L2 


< 


1^1 l~^ + ~ 

X{|?|<i/ 2 }IC|u+(Oexp(-r- 1 ) 


+ 


L2 


X{|CI>i/2}l^|M+(C)exp( 


< ||x{|?|<i/ 2 }IC|exp(-|Cpt)||^cx= I|u+(C)IIl 2 +C>(e *)|||^|u+( 0 I|l 2 

<ori/2)(iK(0iiL^ + iii?i«+(e)iiL^). 


zi±\/iEW. 

2 ’ 


L2 


Here X{.} denote the characteristic functions. The high frequency part decays exponentially and the low 
frequency part becomes the main part. This type decomposition will he used in the proof of the theorem. 
We also remark that the above estimate is optimal without any further assumption for . 

We give a proof of the theorem. First we assume the Cauchy data / = *(wo, mi) is in D{A^) for large 
k and uq £ H. So U{t)f is in as t variable function. We take a smooth cut-off function ip satisfying 


ipit) = 


1, t>2, 
0, t<l. 


We shall estimate the following cut-off propagator as in 0, m, [IS]. 

vit)f = muit)f. 

From (12.21) , we can define its Fourier-Laplace transform by 

1 , 


VfiX) = 


e-^^W{t)fdt 


J-c 

if Im A = —Re iX = —e < 0. Then by Fourier inversion formula, the following identity holds 

1 f 


V{t)f = 


e"^*Vf{X)dX. 


V^. 


ImA= —£ 


By definition, we have the identity 


{dt - A)Vit)f = iP'{f)U{t)f. 

Taking its Fourier-Laplace transform, we get 

{iX-A)V'f{X)=Wf{X). 

if Im A = —£ < 0. Since Al is a dissipative operator the resolvent exists for ImA = —e < 0 and we obtain 

VfiX) = {iX-A)-^P^{X). 

We take its inverse Fourier transform and we get 

Vit)f = ^ / (*A - A)-^A^*^iX)dX. 

\J I'K JlmX— — e 

Since ip’{t) has a compact support, ip'Uf{X) is a holomorphic function and by the results of previous 
section, {iX — A)~^ can be holomorphically extended to Re iX = —a < 0 if a > 0 is a sufficiently small 
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constant and Im iX ^ 0. Since U{t)f is in D[A) x i/, tp'Uf{X) decays sufficiently first as A tend to 
infinity, so we can change the integral contour as follows 


V{t)f = 


/ x(A)(^A-Al)-le*^*V’'t^/(A)dA + ^ 

/ ImA=a V • 


(iA-xl)"^e*^V'CA/(A)dA 


v^. 

= e““‘/i/(<) + hfit) = (High frequency part) + (low frequency part). 


Here x is a cut-off function defined by 


X(A) = 


1, |ReA| > (5, 
0, |ReA| < S, 


for sufficient small S > 0. C which we choose later, is a contour around the origin connecting two points 
iLS + ia . We prove the following bound of Ii by using the method of m which is known as Morawetz’s 
argument 

(4.1) l|e-“‘/i/(t)||«3 < 

First we estimate norm of Iif. Using Plancheral’s identity, we have 

^OO ^ 


/ OO POO -I p 

\\Iif{r)rnodr= ^ e(*^+“)"x(A)(*A - ^)-iV^(A)dA 

-OO J—OO 'y 2.7T 


> —OO 
POC) 


dr 


I-Lq 


< / ||x(A + m)(iA-a-Al) V'H/(A + M)||^^dA. 

J —OO 

By proposition 13.Ill and y cut low frequency, x(A -b ia){i\ — a — A)~"^ are uniformly bounded operators 
on Ho so we estimate 


|/i/(r)||?,„dr<C / ^|J'Uf{X + la) 


' —OO 

^+oo 


^0 


dX 


/ + 00 
-OO 


For last identity, we use Plancheral’s identity again. Since ip' is a compact support function, we have the 
following estimate 


(4.2) 




We give point-wise bound of hf{t) using this estimate. From this inequality, there exist 0 < s < 2 such 
that 

\\hf{s)\\no < CWfWno- 

We define 


dthfit) - Ahf{t) = 


1 




x(A)e(*^+“)‘V^(A)dA + ah fit) 


ImA=a 


= hfit)+ahfit). 

By Duhamel’s principle, we have 

(4.3) hf{t) = U{t-s)hf{s)+ [ U{t-T){l 3 f{T)+ahf{r))dT. 

J S 

We shall estimate the last part. First we obtain 


l|H(t-T)/3/(r)||^^dr = 


(4.4) 


U{t-T) 


^{i\+a)7 


e--'^x{m'Uf{X))dX 


^+oo 


< 


/ ImA=a 

e(*^+“)^x(A)V^(A))dA 


dr 


'Ho 


ImA=a 


dr 


Ho 
























since U{t) is a propagator of a dissipative operator. By Plancheral’s identity we have 

/ +00 ^ 

||x(A + ia)'ip'Uf{X + m)||^^dA 

-OO 

/•+00 

<C U'Uf{X + za)\\ljX. 


We use Plancheral’s identity again and since il)' is a compact support function, we have 

pt p+oo 

I \\U(t-T)hf(TM\l dr <C I 

(4.5) 


Pi ^i-oo 

/ \\U{t-r)hf{r)rp,^dr<C e2“*||V^'(t)C/(t)/||?,„dt 

J S j —OO 

<C^II/ll«o- 

By the Schwartz inequality, from (14.31) . (14.21) and (14.51) . we have 

e"®*||A/(t)||wo < e"®*||C/(t - s)Iif{s)\\Ho + f e-^*U{t - T){hf{T) + aIif(T))dT 

J S 

< C'll/liwo + \\U{t - T)hf{T)\\l^^dT + j^ ||C/(t - T)a/i/(T)||^^dr 

^CII/llH^+C-j /'‘||f/(t-r)/ 3 /(T)||?,„dr+ r ||/i/(r)||^„dr 

L S J —OO 

<C||/||ho- 


1/2 


1/2 


So e “*/i/(t) decays exponentially. 

Next we study l 2 f part. We show the following estimate. 

Proposition 4.3. There exists C > 0 such that for any t > \ and f GTL, the following bound holds 

Whfmno < C{t-^/^f\\n). 

Proof. We choose the contour as follows C = Co UC+ UC_. Co = {Ae*®;s G [—tt, 0]}, C+ = {(1 — s)A + 
s{6 + ia)\s G [0,1]}, C_ = {—(1 — s)A + s{—S + ia)] s G [0,1]}. Here we impose suitable orientation on 
these contours. We estimate each contours in the following integral, 

hfit) = -^ ! (*A - A)-^e^^*f^)fdX. 

VZTT JCoUC+UC_ 

We have ||(iA — A)~^\\n^Ho = on Co by Proposition 13. Ill and the length of Co is 0{\/t). So the 

following bound holds 


'Co 


{iX-A)-^e^^^if'Uf{X)dX 


'Hq 


= o( V/2 sup WWfimn 
E AgCo 


From dUD and ip' is a compact support function, we have 


sup Wip'Uf[X)\\-H = sup 

AgCo AgCo 



e * ^'ip'{s)U{s)fds 


H 


0{l)\W\\M\n 


So we have the estimate for Co part. Next we estimate C+ part, for sufficient large t, we have |A| ^ 
|ReA| ^ s + 1/t on C+. So by Proposition 13.Ill we obtain 

\\{iX-A)-^\\n^no=0{A/^). 


Thus 


'c+ 


{iX-A)-A^^^-tp'Uf{X)dX 


= 0 [ sup \\'ilj'Uf{X)\\n 

Jo agc+ 


In the same way, we can estimate C_ part and we have proved. 


□ 
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Thus we have proved the abstract energy decay of the damped wave equation for sufficiently regular 
initial data. Then by density argument, we obtain the theorem. 

Remark 4.4. In the proof, we essentially use the following facts 

(i) the existence of the resolvent for ImA ^ 0 and ReA > —c, c > 0, 

(ii) uniform estimate of resolvent for large A which comes from |/toA|“^ order estimate of R(X), 

(Hi) order estimate of the resolvent near the origin which comes from |A|“^ estimate of R{X) 

for ReA > 0 and ImA ^ |A| regions. 


5 Proof of Theorem 11.11 

From now on we shall prove Theorem ll.il We recall A = . Since is bounded from 

D{A) to D{A), this operator is also a self-adjoint operator whose domain is D{A). By the spectral 
theorem, we introduce a projection operator f'i^A). Here is a cut-off function defined as follows 


(5.1) 




1, X > e, 
0, X < e. 


By using this function, we define a cut-off operator by 


Then we have 

^B{Af = cfB{A), cfBiAT = 

Since U{t) : H H, we can apply ifsiA) to U{t)f = *{u{t),v{t)), f = ‘(uo,mi) € H hy 

(j)B{A)U{t)f = *{(j)B{A)u{t),(j)B{A)v{t)). 

To prove Theorem 11.31 we use the following decomposition 

U{t)f = cfB{A)U{t)f + (1 - cfB{A))U{t)f 

= (High frequency part) -|- (Low frequency part) 

and estimate each part. First we estimate the High frequency part. 

Proposition 5.1. There exists C > 0 such that for any t > 1 and f € TL, we have the following bound 

\\MA)U{t)f\\n<Ct-^f\\n. 

Proof. Applying similar argument to the previous section, for sufficient regular /, we write 

4>B{A)mU{t)f = MA)-^ [ x(A)(*A - A)-ie*^‘V^(A)dA 

\ JluiX—a 

+ (I)b{A)^ [{iX-A)-^e^^^fyUf{X)d\ 
v27r Jc 

= MA)e-‘^^Iifit) + MA)l2f{t) 

= (High frequency part) -I- (Low frequency part) 

Here 'f’{t) and x(A) are cut-off functions introduced in previous section, a > 0 is a sufficient small constant 
and C is the contour around the origin which is also used in the former section. Since ip{t) = 1 for large 
t, we have 

MA)U{t)f = MA)muit)f 

for large t. So we need the estimate of (j)B{A)if{t)U{t) f. The estimate of e““*/i/(t) is essentially similar 
to the argument to Proposition 14.11 but we need small modification in 631) since we only know the 
estimate (12.2|1 and U{t — t) may be not bounded. In this case, we can prove the following estimate by 
similar argument for arbitrary small e > 0 

f \\e-^^U{t-r)hfrndT <C\\f\\n. 

J S 
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So we conclude the following bound 

(5.2) \\MA)e-‘^*Iifmn = 0(e-(“-^^)‘||/||«). 

Next we shall estimate (/>_b(A)J 2 / by using the operator 

(BA + A)-^ = 

This identity is easily seen from (BA + A) = B^/^(A + From the following identity 

(A^ + BA + A)(BA + A)-^ =Id + \^{B\ + A)-\ 

we also have 


(5.3) B(A) = (BA + A)-i - X^{BX + A)-^R{X). 

Here R{X) = (X^ + BA + A)~^ and (BA + A)~^ are commutative operators. We recall the identity 

(5.4) (A - Al)-i = B(A) = R(X)M(X). 

Applying these identities, we have 

4>B(A)l2f(t) = [ (l)B{A){{BiX + A)-^ + X^{BiX + A)-^R{iX)}M{iX)e^^*^fuJ{X)dX 

V 27 r Jc 

= Jlfit) + J2f(t). 

Here 

Ji/(t) = ^ [ B-i/2^(i)(a +i)-^B-i/2e*^‘M(iA)(/^(A)dA, 

V 27 r Jc 

(5.5) J2f{t) = [ B-i/2^(i)(iA + A)-^B-^/^X^{iX - A)-^e^^*j7uJ{X)dX. 

vStt Jc 

Since 4>{A) is a cut-off operator to the high frequency part. (j){A){iX + A)~^ can be holomorphically 
extended across the origin. Since B^/^ and B~^l^ are continuous on D{A), by interpolation we have the 
following estimate if A is sufficiently near the origin, 

||B ^^'^(j){A){iX + A) = IIa/AB ^^'^(j){A){iX + A) ^u||^ 

= {A(j){A){iX + A)~^u, (j){A){iX + A)~^u)h 

= \\^/Ac|,{A){^X + A)-^urH 

= ll'5^(^)(*^ + ^) 

= 0{\\\fAu\\']j) = 0{\\u\\'\j^ -\- ||u|||f). 

In the same way, using {iX + A)~^A = Id — iX{iX + A)~^, near the origin, we have the following bound 

||B-i/20(i)(iA + A)-^B-^/Hl{iX)\\H^H = 0(1). 


For Jif, we deform the contour to Re iX < 0. Then we have 
1 


\Jim\\H = 


<C f I 

J Im.A=a,|ReA| <.6 


B-i/2^(i)(fA -f A)-^B-^/^A^*M{iX)Jj'Uf{X)dX 


' ImA=a, |ReA| <5 

'"*^‘|||V(^(A)||„dA 


H 


Thus this part decays exponentially. We apply same contour deformation as in Proposition 14.31 to J 2 /. 
Since \\X^{iX + A)~^\\-h^h — 0(A) for Re iX > 0, we obtain \\X^{iX + A)~^\\u^h = 0(l/t) on Co which is 
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defined in the proof of PropositionThen the length of Co is 0[l/t) and (l){A){iX + A) ^ is bounded, 
we obtain 


JCo 


= 0{t 


- 2 | 


WU) 


H 


On C+, we have |A| ~ |ReA| ~ s + 1/t. So ||(iA + A) ^\\-h^h = C){\\\ ^) = 0(t) and 


f B-^/'^(l){A){i\ + A)-^B-^/'^X^{i\ + A)-^e^^*i?UfiX)dX 

Jc+ 

= 0(1) /\(s + lA)2e-“«‘ds||/||„ 

Jo 


H 


= 0{t-^\\f\\n). 

We can get the same estimate for C_. Now by density, we have proved. 

Next we estimate the low frequency part: (1 — (j)B{A))U{t)f = ^B{A)U{t)f. Here we write 


□ 


We assume / is sufficiently regular then by the similar argument used in the previous section, for a cut-off 
function ijj, the following identity holds 

^B{A)i,{t)U{t)f = ^ [ ^B{A)R{iX)MitX)e^’'*^j7uJ{X)dX. 

V^TT JImA—— e 

Writing U(t)f = *{u(t),v(t)), ip'Uf{X) = *{'>p'Uf{X)Q,'ip'U/(X)^), we have the following identity for large 
t 

(5.6) ^B{A)u{t) = ^ / ^B(A)i^(^A)e*^‘ f(^A + B)^f{X\ + V^(A)i) dX. 

v27r JlmA=-£ ^ ' 

We define ipnit) = ip{nt) for the cut-off function ip satisfying ip' ^ 0. Then ip'^ become a Dirac sequence 
i.e. ip'n{t) —>■ S{t), as n —>■ oo where d{t) is the Dirac measure. We change ip in the above integral to a 
sequence of cut-off functions {ipn}- So ip'^Uf{X) —>■ -^^U{0)f = as n —>■ oo and we would like to 

take this limit in the integral. To verify this limiting argument, we use the following lemma. 

Lemma 5.2. If R(X) exists for a X ^ 0 then we have the following bound 

||^B(A)i?(A)u|U = OdAI-^hllff) -f 0(|A|-2(|A| + l)||i?(A)u||ff), 

||^B(A)i?(A)Au||^ = 0(|A|-^||^|U) + 0(|A|-2(|A| + l)||R(A)Au||^), 
for u G H and v G D{A). 

Remark 5.3. If A is bounded, then X^ + XB + A = 0{X^) for large X. So we can expect nearly |A|“^ 
estimate of the high frequency cut-off resolvent. 

Proof. Since (pi^A) : H D{A) and : D{A) — >■ D{A) are bounded, we have 

\\Pb{A)Au\\h = \\^b{A)A\\h^h\\u\\h = \\{^b{A)A)* \\ h^h\\u\\h 
= \\AB^/^iA)B-^^^\\H^H\\u\\H < C\\u\\h. 

So we obtain 

||(<^b(A)(A^ + BX A)u\\h > \X\‘^\\pBiA)u\\H — C'dAI -I- l)||M|l/i-. 

Thus 

Ub{A)u\\h < \X\-^UMA)iX^ +BX + A)u\\h + qA|-^(|A| + l)||w||^^. 

Taking u = R{X)v and u = R{X)Av, we have the desired inequality. □ 
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Since R{\){\ + B) = {1 + R{X)A)/X, by (15.6L we obtain 

=^= [ ^BiA)e^^*-^ilj'Uf{X)QdX 

y ZlT 

+ ^ ( 4i^(^A)AVW7(A)o + R{iX)^{X)^ dX 

y I'K JlmX——£ \'IA J 

= Jl{t) + J 2 (t). 


If supp ijj' C [0, c] we have the following estimate 


IIV'';:^/(A)||« = 


1 

-J= / e-^^^i;'{s)U{s)fds 

V J — oo 


<Ce^'“"|IV''llLi||/ll« 

H 


Ce'='“^|l/||„. 


So for sufficient large t, in the integral Ji, we change the contour to C = C+ U C_ U C. Here C± = 
{±s + is;s G [5, oo)} and C is a contour around the origin. Thanks to the estimate |e*^‘| ^ on 

C+ U C_, if t is sufficiently large, we have 


Ji(t) 


^ [ MA)e^^*\WfWod^ 

\/2tt Jc+uC-UC 


in H. For J 2 , we change the contour to C = C+ U C_ U C. Here C± = {±s + ie; s G [5, 00 )} and C is a 
contour around the origin as in the proof of Proposition |T31 By Lemma l5.21 Lemma 13.51 and Lemma l3.11 
on C±, we have the following estimate, 


-0B(A)i?(zA)HV>'[//(A)o + MA)R{iXWUf[X)^ 


This is integrable and changing the contour, we obtain 


= 0{\X\-^\WUf{X)\\H) = 0{\X\ 


- 2 | 


\n) 


H 


Mt) = 


1 




^^B{A)R{iX)A^f^{X)^ + ^B{A)R{tX)^^{X), ) dX + 0{e 


^ — I 




in H. Thus for sufficient large t, we have 

^BiA)u{t) = ^ ^ (i^(^A)(^A + H)V?t^(A)o + R{iX)^{X),) dX + 0{e-^*\\f\\n). 

Changing ip in the above identity to the sequence {ipn} and we take n to infinity. Since ||'*/'rallLi = 
and supp ^pn C [0,c], the exponential term decays uniformly with respect to n and C is compact, we can 
take limit under the integral sign. So tending n to infinity, we conclude 

^ ^ ^B{A)uit) = ^ f ^B(H)e*^‘(i?(^A)(^A + H)uo+i?(^A)ul)dA + 0(e-^‘||/||«) 

(5.7) 27r Jc 

= Mt)+Oie-^^\\f\\n). 


Finally we apply the identity (15.3|) to get the asymptotic profile. By (15.3F we have 

J3(t) = ^ [ B-^/‘^^{A){iX + A)-^B-^/^e^^*{Buo + ui)dX 
Jc 

+-^ [ B~^^‘^(j>{A){iX + A)~^(iXuo) dX 

2tt Jq 

+ ^ [ B-^/^ip{A){iX + A)-^B-^/^X‘^RiiX)A^*{{iX + B)uo+ui)dX. 

2tt Jq 

In the former integral, we can add contours C± = {±s + is]s G [e,oo)}, modulo exponential decaying 
factor. To the latter integrals, we apply same contour deformation as in Proposition 14.31 Then on the 
contour, (j){A){iX + A)~^ and R{iX) are |A|“^ order. So we obtain these integrals have 1/t bound. Thus 
we get 

Jsit) = ^ [ B-^/^{A){iX + (Huo + ui) dA + 0{l/t){\\uo\\H + ||wi||ff) 

( 5 . 8 ) 27 r Jcuc+uc_ 

= JA^) + C’(l/t)(||uo||/i- + lluillfl-). 
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Since A is self-adjoint operator, it is a generator of an analytic semi-group and we have 


1 

27r 


oUC+uc_ 


(iA + A)-^e^^*dA = e-*^. 


Thus we obtain 


(5.9) 


J 4 (t) (Suo + ui) 

1 f 


27r 


'CoUC+UC_ 


_ ^{A)){iX + {Buo + Ml) dX. 


Since (I5.1|) . 1 — (j){A) = 4){A) cut low frequency, 4){A)(iX -\- A)~^ is analytic near the origin. Thanks to 
the factor, by the contour deformation across the origin, we know the latter decays exponentially. 
By ()5.7|) . (15.811 . (15.911 and density, we conclude 

Proposition 5.4. If </> cut high frequency, then there exists C > 0 such that the following bound holds 
for t > I and u{t) which is the solution of the problem (11.511 with the initial deta *(mo, ui) G B 

11(1 - (j)B{A))u{t) - {Buo + ui) \\h < C(t~^{\\uo\\H + llwillff) + i“^l|wo||ffi)- 


From Proposition 15.11 and 15.41 we have proved Theorem 11.11 

Remark 5.5. From the proof we can get more sharp decay i.e. we can change HvCdMoUff term to 
t~^\\'/Auo\\H- Moreover if we apply the argument of Provosition IHTITI we can get this term is t~°° order. 

Next we show the optimality and an improvement of Theorem 11.11 The following statement says the 
optimality of the decay rate of Theorem 11.11 

Proposition 5.6. Assume that 0 belongs to the spectrum of A. Then 

limsup sup t||M(t) — > 0 

t-^ao ||*('Uo.'tii)ll'H<l 


where u(t) is the solution of the problem (11.511 . 

Theorem 11.11 is an uniform estimate and for individual solution, we can obtain the following result. 
Proposition 5.7. For each solution uft) of the problem (11.511 with *{uo,ui) G %, we have 


lim t\\u{t) - ^ Q_ 


Proof of Proposition^^ From the proof of Theorem 11.11 the main part of u{t) is Jsit). Since {iX + 
A)~^R{iX)iX = 0(|A|) on C, we can apply the same argument as in the proof of Proposition 14.31 
and oi this part is negligible in J 3 . So we have 

I{t) = u{t) - B-^/^e-*^B^/^uo - 

=-^ [ B~^^^^{A){iX + A)~^B~^^‘^e'^^^ (iXuo) dX 

Jc 

■ [ B-^/‘^^{A){iX + A)-^B-^/^X^RiiX)e^^* {Buo + ui)dX + 0{t-^\\f\\n)- 

Jc 


1 

27r 


For latter integral, we apply (15.3F Since R{iX) and {BiX + A)~^ are 0(|A|“^) on C, we can neglect 
the remainder term. For former integral, we apply the similar argument to get heat type asymptotic 
behavior. We obtain 

I{t) = + TT / B-^/'^^{A){iX + A)-^B-^X^{iX + {Buq + ui) dX 


27T 


+ oit-A\f\\H). 

For (j) satisfying dED, we define (t>n{A) = (f>{nA). Then by the spectral theorem, we have 

4>n{A)u ^ u vn H as n — >■ 00 . 
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We decompose I{t) as follows 


+ ^ [ 5-i/20(i)(a + i)-i5-iA2(iA + i)-ie*^‘(?i„(i)5-i/2(5uo + ui)dA 
27r Jc 

+ ^ f 5-i/2<^(i)(iA + i)-i5-iA2(iA + i)-ie*^‘(l-0„(i))5-i/2(5wo + ui)dA 
27 ^ Jc 

+ Oit-^f\\n). 


Since (t)n{A) cut low frequency, (iX + A)(pn{A) is bounded near the origin. So in the first integral, 
the integrand is |A| order and can be estimate Cnt~^ as in the proof of Proposition 14.31 and 15.11 The 
second integral can be estimated as in the estimate of (15.81) and we get uniform estimate. Since 
(1 — 4>n{A))B~'^/'^ (Buo + ui) —>■ 0 in i? as n —>■ iJ. By letting n sufficiently large, this part decays et~^ 
order for any e > 0. Now we recall the fact, tAe~^^u —>■ 0 as t —)> oo, for individual u. This fact is easily 
observed from the spectral theory. Thus we have proved. □ 

Proof of Provosition For simplicity, we can assume uq = 0. From the proof of Proposition 15.71 we 
have 


u{t) - 5-i/2e-‘^5-i/2ui 

1 


27r 


J 5-i/2<^(i)(iA + A)-i5-i/2A25-i/2(iA + A)-^B-^/‘^e^^*uidX + 0{t-‘^\\ui\\H). 


So we get 


{u{t) — B i/^e ‘^5 i/2ui,Mi)_y 


= 7^ /(5-lA2e*^‘(^A + i)-l5-l/2ul,^(i)(^A + i)-l5-l/2ul)HdA + 0(^-2||nl|||,). 
Z-K Jq 


For a € cr(^) Dsupp (f), we take a sequence of functions /„ such that 4’n{A)fn = /„ and ||/n||i? = 1- Here 


tpn{x) 


1 , jx — a| < 1 /n, 
0 , |a: — a| > 1 /n. 


We take u" = B^^^fn, then 

^ = ^\\B-^^^u-\\h < WB-^u^h < C||5-i/2<||^ = C. 

So if necessary, we take a subsequence and we can assume \\B~^Ui\\h converge to Pa > "tjC as n tends 
to infinity. We also note ||A/„||^f is uniformly bounded so from the continuity of 5^/^ and interpolation, 
is uniformly bounded. For A G C, we easily see 

(lA + i)-i5-i/2< - (iA + a)-i5-i/2< ^ 0 


in iJ as n tends to infinity and the convergence is locally uniform. We have 

lim (n”(t) - 5-i/2e-‘^5-i/2n", u^)h = PaJ- [ X^e^^\iX + a)-^dX + 
n^oo ZtT 

In the above integral, we can change C to a simple closed contour C by adding a contour with exponential 
error terms. So we have 

lim {u^{t) - 5-i/2e-*^5-i/2w", u^)h = iPaJ- [ X^e^^*dx{iX + a)-^dX + 0[r'^) 
n^oo 27r 7(5 

= -Pa— [ (2iAe*^* - A2^e*^‘)(^A + a)-ld(^A) + 0(^2) 

2tti Jc 

= /3a(2ae-“‘ - aHe-^A + 

Taking t = \/a and a —>■ 0. Since Pa > l/C, we obtain this integral does not decay faster than order. 
So we have proved Proposition 15.61 □ 
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6 Proof of Theorem 11.3 


Here we study the asymptotic expansion of the propagator for proving Theorem 11.31 In this section, we 
assume H — . We apply the identity (15.31) to (EZD and get 


( 6 . 1 ) 


(j)B{A)u{t) =-^ [ B -\-A) {Buo + ui) dX 

ZTT Jq 

+-^ [ B~^^^^{A){iX + A)~^(iXuo) dX 

27 J' Jc 

+ ^ f B-^^^^{A){iX + A)-^B-^^^R{iX)A^^ {Buo + ui)dX 
27J‘ Jc 

+ ^ [ H-i/20(i)(a + i)-iH-i/2i?(*A)e*^‘(*Auo)dA + O(e-^‘||/||H) 

J c 

= 7 + / + i?i + i?i + Oie-^^WfWu) = I + i + R + Oie-^^WfWu). 


R = Ri + Ri is the reminder term and C is a contour around the origin as in the proof of Proposition 
EM We have estimated I in previous section and get heat type asymptotic behavior modulo exponential 
decay. 

(6.2) I = + 0(e--‘||/||„). 

So we shall estimate I and R. Since (A + H)“^(A + H) = Id, we have 

/= — / B-^/^^{A)B-^/^e^^^uodX- — [ B-'^/^^{A)A{iX + A)-^B-^/^A^*uodX 
27r Jc 27r Jc 

In the former integral, the integrand is holomorphic. So by the contour deformation across the origin, we 
get exponential decay from term. In the latter integral, we can apply similar argument for proving 
heat type asymptotic behavior and obtain the following asymptotics 

(6.3) J = + 0(e-"*||/||«). 

By the assumption, we have the bound 

g-t/2A . 2,1 ^ l 2 ^ C>(t-™). 

Since xe~^, a: > 0 is bounded, tAe~*^ is bounded. So we have 

:L^^L^ = OirA- 

B~^/^ is also a bounded operator on so if uq S L^, using the identity Ae~*^ = , we 

obtain 


(6.4) ||J|U==0(t-”®-i(||no|Ui + ||/||«)). 

Next we study the remainder term R. Since + \ H ^ H = C*(|A|) on the contour, 

using (lOli repeatedly, we obtain, for sufficiently large N 


R=A / e®^*B-i/2^(i)(iA + i)-iH-U2(iA + i)-iB-i/2(Hno + Mi)dA 
27 !" Jc 


+ 


1 

27r 


+ 


1 

27r 



B-^/'^^{A){iX + A)-^B-^X^{iX + A)-^B-^/^ iXuQ dX 

B-^/‘^^[A){iX + A)-^B-^X^{iX + A)-^B-^X^{iX + Ay^R-^^'^ {Buo + m) dX 


+ ...+ [ e®^‘0(|Ar)||/||„dA = + • • • + + 0(r^-i)||/||«. 

J c 

For the latter estimate, we apply same contour deformation as in Proposition 14.31 bv using the estimate 
C>(|A|^) on the contour. Here 


Kn= j + i)-iH-i/2(A2B-i/2(iA + Ay^R-^^AT {Buo + wi) dX, 

Kn=J e^^^B-^/‘^yA){iX + Ay^B-^/‘^{X‘^B-^^‘^{iX + Ay^B-^^AT *Auo dX. 

First we estimate Ki and Ki. 
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Proposition 6.1. There exists a constant C > 0 such that the following estimates hold for any t > 1 
and Mo G n D{\/~A), mi G fl 

\\Ki\\L'i<C{t ™ ^(||Mo||LinL2 + llwilliinis). 

Proof. We use the following identity 

/•i/2 

(6.5) (iA + i)-i = (^A + i)-ie-‘(*A+A)/2 

Jo 


Inserting this identity, we obtain 

Ki=J e^^*/^B-^/^^{A){iX + A)-^B-^X^(iX + {Buq + ui) dX 

f B/2 


dX 


dti e*^(‘-*i).B-i/2^(i)A2(/A + A)-^B-^e-*^^B-^/^ (Buq + ui). 


The first integral is easily treated since we have the (iA + A)~^ = C>(|A|“^) on C. For latter integrals, we 
apply the following identity 

X{iX + A)~^ = —ild + iA{iX + A)~^. 

So we have 

Ki= J e*^*/ 2 C>(l)e-*/ 2 A^-i /2 

r 


-i dX 


dti {Buq + ui) 


f.t/2 


Jc Jo 

r /*^/2 

- dX dh B-^^^^{A)A^{iX + {Buo + ui). 

Jc Jo 

In second and third integrals, integrands are holomorphic with respect to A variable, and t — ti^ tj^^ we 
can deform the contour across the origin and get exponential decay. For the fourth integral, we can get 
heat type semi-group modulo exponential error by applying same argument as in the previous section. 
We obtain 

Ki= J {Buq + ui) 

/•i/2 

+ / i?-l/20(i)i2g-(t-ii)A5-lg-tiA^-l/2 ^ ^ 0(e-"‘)(||uo|U2 + IIM 1 IU 2 ) 

= Ii + h + 0{e ®*)(||uo||l 2 + ||ui||i 2 ). 

Applying same contour deformation as in Proposition 14.31 we get 

||/i||l 2 = (Buo + Ml) IU 2 ) 

= 0(t—l)(||Mo|Ul + ||Ml|UO• 

Since dt\ =0{t) and ||(/)(A)A^e“*/^"^||i 2 _>i ^2 = 0{t ^), the following estimate holds 

II-^ 2 ||l 2 = 0 {t\\ 4 ){A)A^e~^/‘^^\\Li^L 2 sup (Smq + mi) 11 ^ 2 ) 

0<ii<i/2 

= 0{t~^ sup (Bmo + Ml) IIli) 

0<ti<t/2 

= 0{t-^-\\\uo\\m + \\uA\\LA- 

Thus we have estimated ATi part. The estimate of Ki is almost similar. □ 
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We can apply similar idea for the estimate of Kn and Kn terms, as an example, we next prove the 
following estimate 

\\K2\\L^=0{t ^(||wo||LinL2 + llwilliinis)). 

By the identity (I6.5L we have 

K2 = J + A)-^B-^X"^(iX + {Buq + ui) dX 

r /‘i/2 

+ dX dti + A)-^B-y‘^{iX + Ay^B-y^e-^^^B-^/^ {Buo +ui). 

Jc Jo 

Applying the following identity to the second integral 

~ ~ ~ /*i/4 

( 6 . 6 ) {iX + Ay^ = {iX + 

Jo 

we get 

K2 = J e*^‘/2C>(|A|)dAe-*/2^B-i/2 (^Buq + ui) 

r /•i/2 

+ I dX I dh + Ay^B-y‘^{iX + Ayy-*/^^B-y-*^^^B-^/^ (Buo+Ui) 

i/2 pt/4: 


fo 


+ I dX I dh I dta + + 

Jo Jo 


We can apply similar argument as in the proof of Proposition 16.11 and we can get the estimate. 

In the same way, we obtain 

Proposition 6.2. There exists a constant (7 > 0 such that the following estimates holds for any t > 1 
and uo & r\ D{A), ui € LJ C\ 


\Kn\\L'^<C{t "(IliiolliinL^ + llwilliinL^), 

\Kn\\L'^<C{t " "■ ^(lluolliinL^ + ||ui||iinL2). 


By these propositions, we have estimated 4>B{A)ult) part and we have 

(6.7) ~4>B{A)u{t) = B-^/y-^^B-^/yBu^ + u^) + Oir^-yu^B- + hilUi + ||/||„). 

Next we estimate (1 — (j)B{A))u{t) = (j)B{A)u{t) term by sharpening the estimate of Proposition 15.II In 
the proof of Proposition 15.11 we have the decomposition 

MA)mU{t)f = + 0{e-y\f\\u). 

Here tfft) = \ ioi t > 2 and 

yf= f dx f ds B-^/^x^(i>iA){ix + Ayy^^<^*-^^B-^^^y{s)uis)iix-Ayy. 

Jc J — oo 

Since Al is a generator of 17(t), we exchange U{t) to (iA —Al)“^. We apply similar argument to Proposition 
iH By dSSl) and (Ea, we have 

( 6 . 8 ) 

J 2 /= f dX f yis)ds B-^/‘^X‘^(l>iA){iX + Ayy^At-s)B-i/2u(^g'^B-^/‘^{iX + Ay^B-^/‘^M{X)f + R 2 . 

J C J — oo 

Here i ?2 is the reminder term. We recall the support of fj' is in [1,2]. We apply identity (16.51) and J 2 / 
becomes 

(6.9) 

J dxj y{s)ds B-^/'^X'^(l){A){iX + Ayy^At/2-s)^-l/2jj(^^^^-l/2^-^^ ^yl^-tA/2^-l/2j^(^^^j- 
r y /■t/2 

+ dX yis)ds / dh B-^/^X'^yA){iX + Ay y^At-tl-s)^-l/2jj^^^g-l/2^-trA/2 ^-1/2 ^ 

Jc Jl Jo 
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For the first integral, on the contour, we have 


and 

From these estimate, we can estimate the first integral by similar contour deformation as in Proposition 
14.31 and get order decay. In the second integral, we need the following estimate 

||B"i/2e-‘^^/2i?-i/2M(A)/||« = ^ ^)||^||^_ 

This estimate follows from 

+ 1)\\u\\h 

and 

In the above, the first estimate holds by interpolating the following estimate and boundedness on H 

The second one is obtained from the first estimate by taking its adjoint 

= \\{B-^/\-^-^I^B-^/^A^ly\\H^H\\A^/MH 
= WA^/^B-^/^c-^^/^^B-^/^Wh^hW'/AuWh 
= o{t-^/^ + \)\\u\w 

_ -y j'2 

In the second integral, since the integrand is holomorphic and ' is integrable near ti = 0, we obtain 
exponential decay by contour deformation. Using the above estimates, we can repeat similar argument 
for reminder term i? 2 - In this case, one may think the estimate of e~^^'^l'^A is needed but since 

e~tA/'^A = —dter^A!'^ and using integration by parts, we can avoid this term. Since the argument is 
essentially similar, we omit the detail and we get 

Proposition 6.3. There exists a constant C > 0 such that the following estimates holds for any t > 1 
and / G "H 

\\{l-~4>B{A))u{t)\\H<Ct-^-^f\\n- 

From these propositions, we have proved Theorem 11.31 


7 Applications 

As applications of our approach, here we consider some examples of damped wave equation on R"^. We 
study the following divergence form damped wave equation 


(7.1) 



digij{x)dj + a{x)dt 


u{x, t) 


+ P + a{x)dt)u{x, t) = 0. 


Here we write di = and P = — digij{x)dj. We assume gij are smooth functions with bounded 

derivatives and {gijix)) is a family of uniformly elliptic real symmetric matrices, i.e. there exists a 
constant C > 0 such that 

^Id < (gijix)) < Cld 

holds for any x G R"^. Under this assumption P is a positive definite self-adjoint operator on and 
D{\/P) = where and are usual sobolev spaces on R'^. The damping term a (a;) is a smooth 
non-negative function with bounded derivatives. If a{x) > c for a constant c > 0, by Lemma 17.121 below, 
we can apply Theorem 11.11 and Theorem 11.31 Thus we obtain similar decay estimates for the constant 
coefficient case. 
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Theorem 7.1. Let u be a solution of (dD with the initial data u\t=o = uq and dtu\t=Q = ui. We assume 
a{x) > c for c > 0 and define 




)(t) = a{-y 


Uo 


(•)■ 




Ml. 


Then there exists (7 > 0 such that we have the following asymptotic profile 

\\u{t) - v{t)\\L2 < Ct~^{\\uo\\H^ + || mi || l 2 ). 

for Uq € , Ml S and t > 1. Furthermore, if we assume uq G fl and ui G D L^, then we 

also have the following bound for t > 1, 


I|u(t) - vit)U. < Ct-^^/^-yWuoU^nL^ + IlMilU^nii). 


We shall study the case that a{x) may be vanish on a compact set however the geometric control 
condition holds. So we recall about the geometric control condition. For the operator P, we define its 
Hamiltonian flow 4)t{x,£f) = {x{t),^{t)) by the solution of the following the Hamiltonian equation 

{ dx dp d£, dp 

dt d^ ’ dt dx 

(a:(0),C(0)) = (x,0- 


Here the symbol p{x,^) = X]^i=i is associated with P. By the uniform ellipticity, the above 

Hamilton equation has time global solutions for any initial date and we can define the Hamiltonian flow. 

For the damped wave equation, it is known that the mean value of the damping term along a geodesic 
is closely related to the asymptotic profile of the solution c.f. H, m- So we introduce 

1 F 

{a)T{x,0 = —J a{cl)t{x,0)dt 

where T > 0 and we use the notion a(a;,^) := a{x). We study the asymptotic behavior of the damped 
wave equation under the following geometric control condition. 

Geometric control condition; There exist T > 0 and a > 0 such that for all G p“^(l), we have 
{a)T{x,0 > a- 

Under the geometric control condition, we can expect the asymptotic profile is similar to the strictly 
positive case. We first obtain the energy decay. Let U{t) be the propagator of the damped wave equation. 
We have the following energy decay estimate. 

Theorem 7.2. We assume the geometric control condition and a(x) > c for |a::| > R where c > 0 and 
R > 0 are some constants. Let u be a solution of (HID with the initial data u\t=o = uq G and 
dtu\t=o = ui € L^. Then there exists a constant C such that for t > 1 the following bound holds 

d 

I|Vm|||2 + \\dtu\\l2 < C{Y^{gijd,u,dju)L2 + \\dtu\\l2) < Cri(||uo||ffi + IImOIIl^). 


In this case, the damping term is not uniformly positive. But we can prove the similar resolvent 
estimate and we can apply the argument of proving Theorem 4.1. We study the resolvent by dividing 
three range of frequencies. 

For low frequency, we use the following estimate 

Lemma 7.3. If A and 5 > 0 are sufficiently small, then (P + + Aa(-))“^ ■. L'^ ^ Lf exists for ReA > 0 

or |ImA| > (5|ReA| and we have the following bound 

{P + X^+Xa{-))-^ = 0{\X\-y. 


Proof. First we prove the following bound, for sufficient small A > 0 


(7.2) [ iVwpdx + A [ a{x)\u{x)ydx > CX ( \u\^dx. 

JtL‘‘ JRd dRrf 

Taking a smooth compact support non-negative function xix), xix) = 1 for |x| < 1. We set Xrix) = x(rx) 
for r > 0. We have 


||Vu||l 2 > ||XrVM||L2 > ||V(XrW)||L2 - || [V, X^]u|| ^2 . 
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Since a > c for |a;| > R and (Vx)r(a;) = Vx{rx) = 0 for \x\ < 1/r, we obtain 

l|[V,Xr]w||L 2 = r||(Vx)rM||L 2 < Cr ^ a{x)\u{x)\'^dx^ 

for sufficient small r and some C > 0. Since XrW = 0 for |a;| > 1/r, by the Poincare inequality, we have 

llXr^lli^ < C'-||V(XrW)||L2. 
r 

So we get 

||Vu||l 2 > ^r\\xru\\L^ -Cr(^J^ a{x)\u{x)f dx 

Taking r = ev/A for sufficient small e, we have (EH) since a > c for |a;| > R. Now we prove the lemma. 
We use the following identity 

((P + + Aa(-))w,u)L 2 = (P'u,M)i 2 + A^||u |||2 + A / a{x)\u{x)\'^dx, 

From this identity and uniform ellipticity of P, we have 

|((P + A^ + \a{-))u,u)L 2 \> C\\Vu\\l 2 + C|A| / a{x)\u{x)\'^dx - 0(|An|n||i2). 

for Re A > 0 or |ImA| > (5|ReA| if <5 is sufficient small . So for sufficient small |A|, applying (17.211 we get 

|((P + A^ + Xa{-))u, u)i 2 | > C'|A|||m||| 2 . 

From this bound we easily prove the existence of the inverse and the estimate by applying similar argument 
for Lemma ixn □ 

Remark 7.4. The above estimate holds for u G for exterior domain and we can generalize the 

theorem for exterior problem with Dirichlet boundary condition. 

By the geometric control condition, we can control high frequency part and we can prove the following 
statement. 

Lemma 7.5. Under the geometric control condition, P(A) = (P + A^ + Aa(-))“^ : ^ exists for 

X = ir, T G R\{0}, |t| > C for sufficient large C > 0 and we have the following bound 

(P + A2 + Aa(-))-' =0(|A|-i) 

This type lemma is well-known, e.g. H, m, so we omit the proof. From this Lemma and by 
Proposition 13.71 we have 

(A-AI)-i = 0(1 + |A|-1). 

This estimate guarantees the analytic continuation of the resolvent across the imaginary axis for high 
frequency region. 

For middle frequency region, we can use the Fredholm theory. We follow the argument in [TH]. We 
take a smooth compact support function b{x) satisfying, 

c{x) = a{x) + b{x) > c 

for a constant c > 0. Then {P + X^ + Ac(-))“^ exists for pure imaginary X = ir, t G R\{0} and we have 

(P + A^ + Aa(-))(P + A^ + Ac(-))-^ =Id- A6(-)(P + A^ + Ac(-))-^ 

Since (P -b A^ -b Ac(-))“^ takes to and b{x) is a compact support function, by Rellich’s compact 
embedding theorem, b{-){P -b A^ -b Ac(-))“^ is a compact operator. Thus by the Fredholm theory, the 
resolvent exists if and only if the following equation have a solution 

Xb{-){P -b A^ -b Ac(-))“^u = u. 

Taking u = (P -b A^ -b Xc{-))v, the equation can be written as follows 

(P -b A^ -b Xa{x))v = 0. 
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For pure imaginary X = ir, t G R\{0}, multiplying v to the equation, integrating on and taking 
imaginary part, we have 


/ a{x)\v{x)\'^dx = 0. 


So r; = 0 on the support of a. Then by the unique continuation of second order elliptic operators, we 
conclude ?; = 0. Thus for pure imaginary A ^ 0, we have m = 0. So the resolvent exist and by continuity, 
for middle frequency, we have the existence of the resolvent across the imaginary axis. Combining the 
existence of the resolvent for high and middle frequency, we have the following lemma. 


Lemma 7.6. For any e > 0 there exists d > 0 such that R{X) = {P + X^ + Xa{-)) ^ L? ^ L? exists for 
ReA > —(5 if |ImA| > e. 


From these lemmas, we can conclude the energy decay by applying similar argument as in section 4 
c.f. Remark 14.41 


Remark 7.7. We can also prove the energy decay only assuming the geometric control condition, by 
using the argument in see Appendix. 

Remark 7.8. We can also treat a small perturbation of the above case and we can get the energy decay 
for the case that a{x) > 0 does not hold. In the above setting, we take a smooth function b{x) > 0 with 
bounded derivatives and we consider the following damped wave eguation. 

{d^ + P + a{x)dt — eb{x)dt)u(t, x) = 0. 


From the above argument, we get the following bound for X = ir, t £ R\{0} 

||(P + A2 + Aa(-))-'ri||L=<qA|-i||^x|U.. 
Taking u = {P + + Xa{-))v, we get 

|A|||n|U2 <C||(P + A2 + Aa(-))^;|lL2. 


So we have 


||(P + A^ + Aa(-) - eXb{-))v\\L2 > ||(P + A^ + Aa(-))^;||R - C'e|A|||z;|U2 > C|A|||r;|U2 


if £ is sufficiently small. This inequality implies the existence of the resolvent and we can apply previous 
argument to get the energy decay though for applying the previous argument, we need a small modification 
since at first we only know U{t) may be growth like order. 

We also give the asymptotic profile of the solution. We take a smooth compact support function b{x) 
satisfying 

c{x) = a{x) + h{x) > c 

for a constant c > 0. We take 

The following theorem is a generalization of Theorem 11.11 
Theorem 7.9. Let u, v be the solution of the following equations 


(7.3) 

(7.4) 


{df + P + af)df)u = 0, 
u|i=o = uo,dtu\t=o = ui, 

{dt + P)v = 0 , 


on R X R'^ for d > 2, respectively. Under the assumption of Theorem. 17.21 there exists C > 0 such that 
for any uq G , ui G and t > 1, the following asymptotic profile holds 

l|w(0 - c(-)^^^i'(0IIl2 < Ct~^{\\uo\\m + I|mi||l 2 ). 
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To prove the theorem, we need the asymptotic expansion of the resolvent near the origin. Since c(x) 
is uniformly positive, the resolvent R{X) = {P + + Ac(-))“^ exists near the imaginary axis. We have 

the following identity 

(7.5) (P + + \c{-))-\P + A^ + Aa(-)) = Id + X{P + A^ + \c{-))-^b{-). 

We use the following estimate for the cut-off resolvent. 

Lemma 7.10. Let d > 2 and x is a continuous compact support function then there exists a constant 
5 > 0 such taht the following estimate holds for |A| with ReA > 0 or |ImA| > 5|ReA|, 

\\xRWx\\L^^L^=oii). 


Proof. By uniform ellipticity, we have 


\{{P -f A^ -f Xc{-))u,u)l2\ > CIIVmIIIs -f CIA] [ c{x)\u{x)\‘^dx - C\X\‘^\\u\\‘l2 
for ReA > 0 or |ImA| > (5|ReA| if S is sufficient small. Since c{x) > c, for sufficient small A, we have 

|A| / c{x)\u{x)\^dx-mu\\l 2 > 0 . 

If d > 2, by the Hardy’s inequality, we obatain 

Here (x) = vT+~?. So we have 

||((a;)(P + A^ -f Ac(a;))'u||i2||(a;)"^M||z,2 > \{{x){P -f A^ -f Xc{x))u, {x)~'^u)l2\ 

= \{P + X'^ + Xc{x))u,u)l^\ > C'||(x)“^m| 1 | 2 . 

Taking u = R{X){x)~^v, we obtain 


(7.6) 


{x)-^R{X){r^v\\^2<-\\vh2. 


From this inequality, we have the estimate 

For d = 1, 2, we can give the following results. 


□ 


Lemma 7.11. If x is a real-valued continuous compact support function, then there exists a constant 
d > 0 such that the following estimates hold for sufficient small |A| with ReA > 0 or |ImA| > (5|ReA|, 


||xP(A)xl|L^^L^= 0 (log|A|) 

\\xRWx\\l^^l^=o{\x\-^/^) 


(d = 2), 
(d=l). 


Proof. By uniform ellipticity, we have 

(7.7) |((P + A^ + Xc{-))u,u)l^\ > C{\\Vu\\l2 

We prove the following estimate. 


klli2)>q|(|v| + |Ar/>|| 


(7.8) 


ii(ivi + iAr/>iii2>{ 

I C\X\^'qxu\\i2, (d=l). 


2 

L 2 - 


We take a compact support cut off function x(^) which satisfies x(^) = 1 if |5| is sufficient small. Then 
by the Plancherel theorem and (|^| -f |A|^/^)“^(1 — x(?)) bounded for sufficient small A, we have 

||x(|V| + |A|V2)-1«|1|, < ||;^(|V| + \X\^/T^x{D)u\\l 2 + ||(|V| + |A|V2)(1 _ x{D))u\\l 2 

< ii(ivi + \x\^/^)-^xiD)u\\i^ + iKiei + iAr/^)-i(i - m)mh 
<m\ + \M^^^)-^mMh+c\\u\\i2 
<ll(iei + |Ar/T'x(C)lli2||«lli^ + qHli- 
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By easy computations, we obtain 


Taking u = (|V| + |A|^/^)u, we have (Tra . From (I7TD and (Tra . we get 

Clog|A||((P + + Xc{x))u,u)l 2\, {d = 2), 

C\X\~^/^\{{P + A^ + Xc{x))u,u)l 2\, {d = 1). 

Now we take u = (P + A^ + Xc{x))~^x'v then we have 

(d = 2), 

(d=l). 

Dividing by ||xu||l 2 , we have proved. □ 

We take the inverse of the right hand side of the identity (Ha. We define 

P(A) = VK^X{P + A^ + Ac(-))-'\AO- 

Since b(x) is compact support function, by the above lemmas, we can assume ||P(A)||i 2 _>i 2 < 1 for 
sufficient small A. For such A, we can construct the inverse by Neumann series. 

{Id - X{P + X^ + Xc{-))-Hi-))-^ 

CO 

= Id+ ^(-1)"A(P + A^ + Ac(-))-'xA( 0{AP(A)}" v^. 

n—0 

Thus we have the following identity for the resolvent P(A) = (P + A^ + Xa{-))~^ and Rc{X) = (P + A^ + 

Ac(-))-i 

P(A) = (P + A2 + Ac(-))-^ 

OO 

(7.9) + ^(-1)"A(P + A^ + Ac(-)) ^y/b{-){XR{X)}"‘\/b{-){P + A^ + Ac(-)) ^ 

n—0 

= Rc{X) + XRc{X)b{-) Rc{X) + Ri{X). 


II II II \2 _L w tt-i II C'logp' ^ C'log P' UH XuLl2, 

XM + A + Ac(a:)) Xv\\l^ < < a ... ii 

[ C A ^'^\{{v,xu)l^ <C\X\ ^'^\\v\\lA\xu\\l^, 



Here Pi (A) is the reminder term. So we have the asymptotic expansion of the resolvent. We shall prove 
Theorem 17.31 If d > 2 then P(A) is bounded so we have ||AP(A)||i 2 _>i ^2 = 0(|A|) and ||Pi(A)||i 2 _>i 2 = 
0(1) near the origin. To the first term Pc(A), we can apply similar argument in Theorem 11.11 and get 
the heat type asymptotic. Pi (A) term becomes lower order term since ||Pi(A)|| = 0(1) is lower than the 
main order |A|“^ and we can get order decay. The main difference is the second term and we study 
the term. We have ||())c(P)APc(A)&(a:)Pc(A)|| = 0(1) if (j)c{P) cut low frequency so, corresponding to 
Proposition 15.11 we have order decay to the high frequency part. For low frequency part, we follow 
the proof of Theorem ll.il To obtain (15.71) . we do not use the identity (15.31) and we only use these estimate 
of the resolvent which is similar in this situation. Thus the main term is J 3 and we have 

(7.10) u{t) = ^ ( MP)e^^* {R{iX){iX + c{-))uo + R{iX)ui) dX + 0(t-p. 

Jc 

Inserting dzH), since Pi (A) and the first term are a lower order, we get 


(7.11) 


u{t) = ^ [ ^c(P)e*^* (Pc(iA) + iXRc{iX)b{-)RciiX)){c{-)uo + ui) dX + 0{t-^) 
Jc 

= Li + ^P2+0(rp. 


In the above integral the first term is similar to the strictly positive case, so we can get the following 
asymptotic profile by applying similar argument to prove Theorem ll.il 


Pi = c(-) *^c(-) ^/^(c(-)uo + Ml) + C>(^ ^)- 


To estimate P 2 part, we recall the following properties of our heat propagator. 
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Lemma 7.12. There exists a constant C > 0 such that the following bounds hold for u G D and 
t > 1 


\\e-*^uh^<C\\u\\L^, 


Proof. Probably this result is well know c.f. [5] and we give an outline of the proof. We use the following 
Nash’s inequality 

Since P is uniformly elliptic and c(-) is strictly positive bounded function, this inequality is easily obtained 
from the following usual Nash’s inequality 


Let u{t) = u{t, ■) = e then by the integration by parts, we have 

■ / c{xY^‘^u{t,x)dx = — / di(gij{x)djc{x)~^/^u{t,x'^dx = 0. 


_d 

dt 


Thus the propagator preserves this modified total heat. Under our assumption, also preserve 
positivity see [S]. So we can only consider the positive solution and we assume u(t) is positive. From the 
conservation of the modified total heat, we have the bound. Next we define 


H[t) := 


Wuimh 

H-y/MmV 


Since the modified total heat is preserved, by the Nash’s inequality, we have 



2(Pu, u) 




This inequality implies 
So we obtain 


-H{t)-'^/‘^ < < -Ct. 


\u{t)\\L^ < Ct-'^/^WuoW 


Li- 


From this estimate, we can prove the estimates by using duality argument. For positive v and s < t, 
we take 

v(s,x) = v(s) := 


Then we have 
d 
ds 


j u{s,x)v{s,x)dx = j (^—Pu{s,x)^v{s,x)dx + j u{s,x) (^Pv{s,x)^ dx = 0. 
Thus the above integral does not depend s so taking s = t and s = t/2, we have 


/R'i 


u{t,x)vdx = / u{tl2,x)v{tl2,x)dx < \\u{t/2)\\L2\\v{t/2)\\L2 < Ct '^^‘^\\u\\l 2 \\v\\l^ 




By this inequality and duality, we get the L°°Lf estimate. 

Using this bound, we estimate L 2 in the decomposition (17.101) . From the identity (15.31) . we have 
(7.12) i?c(A) = c(-)”^/^(A + P)-^c{-)-^/‘^ - A2c(-)”^/^(A + P)-^c{-)-^/‘^Rp\). 

Inserting this identity to L 2 , since A^c(-)“^/^(A + P)-ic(-)“^/^i?c(A) = (!1(|A|), we obtain 


□ 


U2 = y (^c(U)e*^UAc(-) ^/^(A + P) ^c(-) ^/^(A + P) ^/^(c(-)uo + rti)dA + C>(t 

= L^ + 0[t-^). 
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So we estimate L 3 part. We use identities (16.51) and (16.61) . By these identities, we have 

L 3 = y + P)-ic(-)-^^"(c(-)uo + Mi)dA 

r rt/'i 

+ / ^c(-P)e*^*tAc-i/2 / e-‘i(*^+^)dtic-i/26(.)c(-)-^/"(iA + P)-ie-‘(*^+^)/Vi/2(c(.)uo+ui)dA 

Jc Jo 

r ~ /**/4 

Jc Jo 

f /**/2 - /*i/4 

+ / 0 e(P)e*^‘^Ac(•)■^/" / + ui)dA 

Jc Jo Jo 

= _L4 + + iv 7 

Since 6 (a:) is compact support function and d > 2, from Lemma l7.12[ we have 

||g-tP/2c(2.)-l/2^(2.)c(2;)-l/2g-tP/4||^2_,.^2 < \\e~*^^'^\\L^^L^\\c{x)~^^'^b{x)c{x)~'^^'^ 

< 

From this estimate, we know L 4 is a lower order term. For L 5 , since 

f.t/2 




= max{l,t"'^/‘‘+^+^} 


pi-s-ps 


for any e > 0 , we have 

t/2 


< 


< 


10 


P^(3‘-4‘i)/4(ac-i/2e-‘iPc-i/26(.)c(-)-i/"(*A + P)-ie-‘P/4c-i/2dti 

L2-S.L2 

c-i/26(-)c(-)-'/"e-‘^/4*A(zA + P)-ic-i /2 


giA(3t-4ii)/4^-l/2g-tiP 


Li-i.L2 


L2_>li 


dti 


'0 


gZA(3t-4ti)/4^-l/2g-tiP 


Li->.L2 


dti 


c ^/^ 6 (-)c(-) ^/^lliille *^/^||L 2 ^icx,||iA(iA + P) 


L2_>L2 


< C'max{e-3‘i“V^e-““^/^}max{l,t-‘^/'‘+i+"}t-'^/^ 


From this estimate the integrand is 0(1) as A variable and by contour deformation as in the proof of 
Proposition 14.31 we get further decay. Thus we have 

IITsIIl^ < Cmax{l,t-'^/4+i+4r'^/4-i < Ct-\ 

We can apply similar argumetn to estimate Lq term and we get order decay. For P 7 term, the 
integrand is holomorphic so thanks to the factor, we can get the exponential decay estimate, 

by contour deformation across the origin. Thus we have proved Theorem 17.91 For d = 1,2, in the same 
way, we can also prove the following result. 

Theorem 7.13. Under the assumption of Theorem \ 7 .21 for any e > 0 and d = 1,2, there exists C > 0 
such that for t > 2 the following estimates between the solution u of the equation (1131) and the solution 
V of the equation (El hold 


l|w(t) - c(-)^/^p(t)||p 2 < Ct ^ logt(||uo||pi + ||ui||i, 2 ) 
\\uit) - c(-)^/^p(t)||L2 < Cf^/^dluolliJi + ||ui||i2). 


(d = 2), 
(d=l). 


Remark 7.14. In the above theorems we have the arbitrariness to the choice of c(-) term. Probably our 
estimates come from the arbitrariness of c{-). For simplicity, we consider the following heat equation 

((1 + a{-))dt - A)u = 0, u|t=o = (1 + a{-))~^uo 

Here a(-) is a non-negative compact support function. Then by Duhamel’s principle we have 


i{t) = e‘^u( 0 ) - [ e^*-^^^a{-)dsu{s)ds. 

Jo 
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By integration by parts, we obtain 


u{t) = ^^^a{-)dsu{s)ds — j (^dte^* a{-)u{s)ds — e^^‘^^a{-)u{t/2). 

In the right hand side, the first term is the top term. Since a{x) is a compact support function, we have 

< ||e‘/^^|Ui^i2||a(.)||Li||4i/2)||L^^L== < Ct-^^/^umiL^- 
From well-known estimates of heat propagators, we get the following estimate 


i/2 


g(i ^l^a{-)dsu{s)ds 


L 2 




t ^||uo||l 2 (d = 2 ), 

(d=l). 


We can treat the another integral by same manner and get similar estimate. So the above estimate seems 
almost optimal. 

In the same manner, we can generalize Theorem ll.3l 

Theorem 7.15. Let u, v are the solution of the equations (1731) . (1731) respectively. Then under the 
assumption of Theorem 7.1, for any e > 0, there exists C > 0 such that for any uq € D and 
{x)uo € Lfi, ui G n and t > 2, the following asymptotic profile hold 


\Ht)-c{-)^/MmL^<ct 


-l-d/4 


(l|wo||ifinLi + llwilli^nLi) 


Proof. The proof is essentially same and we give an outline of the proof. We use the expansion (1731) . Rc 
term can be treated by similar manner in the proof of Theorem 1 1.31 and we get the heat type asymptotic. 

□ 


8 Appendix 

In this appendix, we prove the energy decay estimate under the geometric control condition. We estimate 
the resolvent dividing three frequency parts in the same way as Theorem [721 We can use Lemma 17751 to 
estimate the high frequency part. So we need to estimate middle and low frequency parts. We take 

L = P + A^ + Aa(-) 

and estimate ||Lm||j ;,2 from below. The following argument is essentially similar to [7]. By the geometric 
control condition, we can take a sequence (xn) C such that a(xn) > a/T and any point of is 
at bounded distance from the set U„{a;„}. By the uniformly continuity, we can find r,/3 > 0 such that 
a(x) > fi > 0 on cj = U„i?(a;„,r). 

We estimate the low frequency part. First we assume |ReA| < |ImA|. By taking imaginary part of 
{Lu,u)]^ 2 , we easily have 

(8.1) II^m||^ 2 + e^|A|^||u||^2 + C'e|A|^||M||^2 > -®^||v^u||^2. 

We can take e > 0 arbitrary small and C > 0 is some constant. As in [7], we set 

+A2)e-W'‘ 

here h > 0 is sufficiently small. Then we can take a uniformly bounded weight tp such that the following 
sub-elliptic estimate holds 

(8.2) WQhuWl. > Ch\\u\\l. 

here C > 0 and u satisfies = 0, see [7] for the proof. We take a real-valued cut-off function y G 
such that x(s) = 0 for s > (3 and y = 1 in a neighborhood of 0. Then = y o a vanishes on u and any 
derivatives of is controlled by a. We set v = and compute 

QhV = h^e^^^a'^Lu — hfe^^^Xa^^au -\- hfie'^^’^ [P, a^] u. 
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Since v\ai = 0, from this identity and (18.21) . we obtain 

\\h^e^/^a^Lu\\l 2 + Ilh^e^/'^Aa^onlli. + [P,a^]u\\l 2 > WQhvWh > Ch\\v\\l.. 

Since (p is uniformly bounded, for fixed /i, is uniformly positive and bounded . Thus we have 

\\Lu\\l. + CmauWl. + C|| [P, a'^] u\\h > ^||a%||i.. 

Since derivatives of a° are controlled by a, for sufficient small |A|, we have 

(8.3) e^/^\X\‘^\\Lu\\l2 + Ce^/^\X\‘^\\^u\\l2 + ^ 

*j=i 


For fc = 1,2, • • • ,d, using integration by parts, we have 

d d 

{Pu,{dka'^)^u)L2 = ^ {gij{dka'')diU, {dka‘')dju)L2 + 2 ^ {{dka‘')diU, g,j{djdka‘')u)^2. 
i,i=l 

Since derivatives of a'^ are controlled by a, from this identity, we have 

d 

\\Lu\\l2 + \\y/au\\l2 > C* XI \\idka‘')^iu\\l2 

fc,2— 1 


for sufficiently small |A|. So from (j8.3p . we get 

(8.4) e^/^|Ap||Lu|||2 + C'e^/^|Ap||-\/au|||2 > —^ ^ ||a°u|||2. 

Thus from (18.ip and (18.4p . we obtain 


||Lm||^ 2 + C'e^|Ap||u||^2 + C'e^/^|Ap||-v/au||^2 + C'e|Ap||M||^2 > ll^/aw||^2 


e3/2|A|2 

C 


||a°u|| 


2 

L 2 - 


Now taking e sufficiently small and |A| < e, since a°(x) = 1 if a{x) = 0, we conclude 


(8.5) 


\\Lu\\l2>Cmu\\l2 


From this estimate, we get P(A) = 0(|A| ^) for sufficiently small |A| around the imaginary axis. If 
ReA > |ImA|, by taking real part of {Lu,u)]^ 2 , we have 

||Lm||^ 2 + £^|Ap||u||^2 > Ce\X\^\\^/au\\\2■ 


By similar argument, we can also prove (El) in this case. So we have (18.51) and get P(A) = OdA] ^) for 
ReA > 0. Thus we have estimated the low frequency part. The estimate of the middle frequency part 
is essentially same and more easy (we only need the reversibility of P — A^ + iXa if A is real and in a 
compact subset of (0,oo)) so we omit the proof. From these estimate, We conclude the following energy 
decay estimate. 

Theorem 8.1. We assume the geometric control condition. Let u be a solution of El with the initial 
data u\t=o = uq G and 9tM|t=o = ui G L'^. Then there exists a constant C such that for t > 1 the 
following bound holds 


< C'( X 9ju)i2 + ||9tu||^2) < Ct ^(llwollffi + ||mi)|| 
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